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We report on the application of a dynamic decoherence control pulse sequence on a nuclear 
quadrupole transition in Pr^^ :Y2Si05 • Process tomography is used to analyse the effect of the 
pulse sequence. The pulse sequence was found to increase the decoherence time of the transition to 
over 30 seconds. Although the decoherence time was significantly increased, the population terms 
were found to rapidly decay on the application of the pulse sequence. The increase of this decay rate 
is attributed to inhomogeneity in the ensemble. Methods to circumvent this limit are discussed. 
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PACS numbers: 03.67.Pp,76.70.Hb,76.30.Kg,76.60.-k,76.70.Hb,76.60.Lz,76.60.Es 



Recent experiments have demonstrated that optically 
active centers in solids and in particular rare-earth ions 
in crystals show promise as an alternative to atomic 
based ensembles for quantum information processing 
applications P, I^. Rare-earth doped ions provide multi- 
level systems possessing optical transitions with long co- 
herence times and in crystalline hosts, high optical den- 
sities are achievable. These solid state centers have the 
significant advantage over atomic systems that they are 
stationary, thus allowing spatial and frequency variation 
in an ensemble induced by interactions with light fields to 
be maintained for appreciable times. This last attribute 
has already been exploited extensively in the field of time 
domain optical processing, where complex temporal and 
spatial relationships between different interacting light 
fields are employed to process classical information, in a 
manner not possible in atomic systems 3j . 

A common feature of many of the proposed ensem- 
ble based quantum information processing applications, 
including quantum memories, single photon sources and 
quantum repeaters, involves the storage of quantum in- 
formation on long lived low lying states 0, y, 0|. In the 
case of rare-earth centers it is proposed to store the quan- 
tum information on the centers' ground state hyperfine 
transitionsI3,0. The duration and fidelity that such in- 
formation can be stored depends on the coherence time of 
the hyperfine transitions and, therefore, it is important to 
maximize these coherence times. The coherence times of 
hyperfine transitions in rare-earth centers are typically of 
the order of ms but in recent work we have demonstrated 
ways where these can be extended and have achieved a co- 
herence time of 82ms in the case of mj = —1/2 ^ +3/2 
transition in In this letter we demon- 

strate a further dramatic increase of the coherence time 
to several tens of seconds through the use of a Bang Bang 
dynamic decoupling sequence. 

Dynamic decoupling methods for open systems were 
initially developed for use in NMR spectroscopy to selec- 
tively remove contributions to the spin Hamiltonian [Tol 
ITU . There is growing interest in applying these tech- 



niques and in particular ''Bang Bang^' decoupling se- 
quences to quantum information systems to decrease 
their rate of decoherence 12, 

3 El E Ei. "Bang 
Bang" sequences operate by decoupling the quantum sys- 
tem from the bath through the application of a periodic 
control Hamiltonian. In the current work this is imple- 
mented by a series of pairs of hard pulses (tt, — tt) sep- 
arated by the cycling time Tc. The Bang Bang pulse 
sequence can theoretically rephase all the coherence in 
the quantum system, thereby making T2 = Ti if the fol- 
lowing criterion |12j is met. 
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where uOc is the cutoff frequency of the dephasing bath 
and CTc is the period between the tt pulses in the Bang 
Bang sequence. The pulses are assumed to be 'hard' such 
that during the pulse any evolution of the state other 
than the action of the driving field can be assumed to be 
negligible. 

For Pr ions substituting for Y in ¥28105 the main 
source of decoherence of the its ground state hyperfine 
transitions is magnetic interactions with nuclei in the 
host0, 0. Y possess a nuclear spin and has a mag- 
netic moment of jy=209Hz/G. The other nuclear spin 
in the host is '^^Si found in natural abundance (~4%) 
with a moment of 84511 z/G. The magnetic field seen 
by any given Pr ion fluctuates over time due to reso- 
nant cross relaxations between the host spins. The sen- 
sitivity can be calculated using the reduced Hamilto- 
nian of the electronic ground state hyperfine structure El- 
Praseodymium has a nuclear spin of 5/2 and as a result 
of a pseudoquadrupole interaction its electronic ground 
state splits into three doubly levels with splittings of the 
order of 10 MHz. The final degeneracy can be lifted by a 
magnetic field through an enhance nuclear Zeeman inter- 
action. In the region where the applied field produces an- 
ticrossings a field direction and magnitude can be found 
such that there is no first order Zeeman shift for the 
mj = —1/2 ^ +3/2 hyperfine transition, demonstrated 
in recent work0. It has been shown that the dephasing 
perturbations acting on the Pr ions are predominantly 
occurring on a timescale longer than ^lOms 9j. Given 
that Rabi frequencies greater than lOOkHz are readily 
achievable the period between the applied pulses in the 
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FIG. 1: Pulse sequences used in experiment: a) 2 pulse spin 
echo, b) Inversion Recovery Ti c) Bang Bang pulse sequence. 
The coherence generated by the initial pulse is rephased by a 
pair of TT, — TT pulses separated by the cycling time Tc, iterated 
N times. The laser is common to all pulse sequences 




FIG. 2: All ground state hyperfine levels other than mi = 
— 1/2 interact via RF frequencies lji-3 with laser radiation ei- 
ther at the read frequency cOr or pump frequency oOp. Through 
spontaneous emission from the electronic excited state the 
group of ions interacting with both Ur and Up via any com- 
mon excited state will be holeburned into the mi — —1/2 
state. Bold states are the Critical Point transition. 



Bang Bang sequence can be made such that the inequal- 
ity (1) is satisfied for the frequencies of the dominant 
perturbations. 

As in previous work a Pr^^ '.¥28103 0.05% concen- 
tration crystal maintained at temperature of ^1.5K in 
a liquid helium bath cryostat 9j. Yttrium orthosilicate 
(¥28105 ) is a low symmetry host with two crystallo- 
graphically inequivalent sites where Pr can substitute for 
Y, labelled 'site 1' and 'site Site 1 ions are used in 

this work. The coordinate system chosen is the C2 axis 
is z is the direction of the predominate polarization of 
the optical ^H^ — transition and x is perpendicular 
to both. 

The magnetic fields to achieve the Critical Point field 
configuration (^cp = (732, 173, -219) G 9]) were sup- 
plied by two orthogonal superconducting magnets sup- 
plying a z field, and x, y field. The sample was rotated 
about the z axis to provide the correct ratio of fields along 
the X and y axes for the critical point in magnetic field 
space. The field in the x, y plane could also be adjusted 
using a small correction coil mounted orthogonal to the 
main x^y. The inhomogeneity in magnetic field across 
the sample was measured using a hall probe to be <2G. 

Raman heterodyne was employed to investigate the 



ground state hyperfine transitions using an experimen- 
tal setup similar to that described in previous work0. 
The experiment was performed using a Coherent 699 fre- 
quency stabilized tunable dye laser tuned to the ^H^ — ^D2 
transition at 605.977nm. The laser's frequency was sta- 
bilized to a sub kilohertz linewidth. The laser power inci- 
dent on the crystal was 40m W, focused to ^100 jam and 
could be gated using a lOOMHz acousto-optic modula- 
tor. The hyperfine transition was excited using a six turn 
coil with a diameter of 5mm, driven by a 101^ RF ampli- 
fier resulting in a Rabi frequency Vtrf = 91kHz. The RF 
pulse and digital control sequences were generated using 
a direct digital synthesis system referenced to an oven 
controlled crystal oscillator. The pulse sequences used in 
the experiment are illustrated in Fig. ^The Raman het- 
erodyne signal, seen as a beat on the optical beam, was 
detected by a 125MHz photodiode. This signal was ana- 
lyzed using a mixer and a phase controlled local oscillator 
referenced to the RF driving field. 

At the Critical Point field the m/ = —1/2 ^ +3/2 
transition was observed at 8.646 MHz with a inhomoge- 
neous linewidth of AkHz. This linewidth was found to be 
insensitive to changes in the magnetic field of the order 
of -lOG. 

Prior to applying each Raman heterodyne pulse se- 
quence the sample was prepared with the optical/RF 
repump scheme as shown in Fig. [21 The repump fre- 
quencies were uor — = l8.2MHz^ uji = 12.2, uj2 = 
15.35MHz and cjs = 16.3MHz. The repump RF was 
pulsed with a duty cycle of 10% to reduced sample heat- 
ing, while the repump laser frequency uOr was scanned 
200kHz to hole burn a trench in the inhomogeneous op- 
tical line where detection would take place. This repump 
scheme ensures that all Pr ions interacting with the laser 
radiation are forced into the m/ = —1/2 state, creat- 
ing a pure state ensemble. It also ensures there is no 
initial population near the laser frequency used for Ra- 
man heterodyne detection. The use of a sub kilohertz 
linewidth laser and the repump scheme resulted in a sig- 
nificant improvement in the signal to noise compared to 
work performed earlier [^J. 

Figure shows the decay of the amplitude of two pulse 
spin echoes as a function of the delay between the pulses. 
Two data sets are shown, the first is for an applied mag- 
netic field optimized for long decay times and the second 
for a field detuned by 5G in the z direction from the op- 
timal value. The significantly longer decay time for the 
optimised field in the present work, compared to previous 
work 9j is attributed to better identification of the crit- 
ical point enabled by the addition of the correction coil 
enabling significantly more precise magnetic field adjust- 
ments. Besides being longer the decay can no longer be 
described by standard echo decay function with a sin- 
gle time constant |I9| . There are three distinct regions. 
For pulse separations less than 20m5 the decay rate is 
less than 1/4^"^. At 30ms there is a distinct shoulder 
with the decay rate increasing to 1/0. 4s~^ as the pulse 
seperation reaches 60ms. From 150 ms onwards the de- 
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FIG. 3: Spin echo decay at the Critical Point field (o), and 
detuned from the Critical Point field by ^2G (•) and 5G (o) 
in the Z direction) 



cay rate asymptotically approaches a value of l/0.86s~^. 
This asymptotic decrease in the decay rate was only ob- 
served for magnetic fields within 0.5G of the optimal field. 
When the field was more than 0.5G away from the crit- 
ical point a simple exponential decay was observed for 
delays longer than 50 ms. 

The shoulder in the decay at 30ms is interpreted as 
indicating that the majority of the dephasing is due to 
perturbations that occur on time scales between 10 and 
100ms. The asymptotic behavior of the decay is at- 
tributed to a variation in the T2 within the ensemble 
resulting from inhomogeneity in the magnetic field across 
the sample. Ions experiencing a field closer to the Critical 
Point condition will have a longer T2 and consequently 
their contribution to the echo intensity will dominate for 
long pulse separations. The asymptotic decay rate there- 
fore gives an upper limit for the contribution to the deco- 
herence rate due to second order magnetic interactions. 

The Bang Bang pulse sequence was investigated us- 
ing an initial delay of ri = 1.2ms and varying the cy- 
cling time Tc from 20ms to 0.5ms, as shown in Fig. ^ 
Also shown in Fig. ^s the result of an inversion recovery 
measurement used to determine the lifetime of the tran- 
sition. The inversion recovery measurements were per- 
formed using the pulse sequence described in Fig. ^b). 
As can be seen in Fig. ^the Bang Bang sequence sig- 
nificantly increases T2, though for the shortest value of 
Tc T2/T1 < 1/4. In Figure ISl the coherence times T2 for 
each of the data sets are plotted as a function of the cy- 
cling time Tc . This shows that for Tc < 5ms ^ there are 
significant gains in T2 made by reducing Tc, while reduc- 
ing Tc further only slightly increases in T2. Also shown 
in Figure [51 is the same measurements made with the dc 
magnetic field detuned so as to increase the transitions 
magnetic field sensitivity. The field was detuned such 
that the two pulse echo T2 was reduced to 100ms. Al- 
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FIG. 4: Bang Bang decoupled echo decays with Tc = 
7.5ms 10ms (o), 15ms (V), 20ms(x) corresponding to T2 = 
27.9,21.1,15.2,10.9s. Inversion recovery measurements (o) 
yield Ti = 145s 
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FIG. 5: Dependence of decoherence time on the Bang Bang 
cycling time Tc both at the critical point (o) and with the 
magnetic field misaligned to give a coherence time of T2 = 
100ms (x). Trend lines are included but do not represent a 
physical model. 



though for long Tq there is a significant reduction in the 
T2 observed using the Bang Bang for the measurements 
made with the detuned field, T2 in the limit of short Tc 
appears to be the same for the two magnetic field condi- 
tions. This suggest that the residual decoherence in the 
limit of short Tc is not predominately due to magnetic 
field fluctuations. The residual decoherence is possibly 
due to detuning and pulse area errors in the decouplng 
sequence. 

To assess how well arbitrary quantum state are pre- 
served by the Bang Bang pulse sequence we performed 
process tomography on the input state and for 1, 10, 
100 and 1000 iterations of the Bang Bang cycle. The 
initial delay was ri = 1.2m5, with a cycling time of 
tauc = 2ms. The total period over which the tomog- 
raphy was performed was ~4ms, ~40m5, ~400m5 and 
~45 respectively. The imaginary component of the pro- 
cess tomography is only shown for the input state since 
it never contributed more than 10% of the signal. Ideally 
the Bang Bang process operator is the identity matrix, 
leaving the state unchanged. It was observed that the 
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FIG. 6: Process tomography of the Bang Bang pulse sequence 
for 1, 10, 100 and 1000 iterations for both experiment and 
theoretical modeling. The imaginary component of both the 
experiment and modeling is always zero and is omitted for 
clarity 



component of the Bloch vector in the coherence plane 
for a given state was preserved well, while the population 
component of the Bloch vector rapidly decayed. The fi- 
delity of the Bang Bang process for 1, 10, 100 and 1000 
iterations was 99%, 65%, 54% and 43% respectively. 

The evolution of the ensemble was modeled using the 
Bloch equations assuming an infinite Ti and T2, with an 
inhomogeneous linewidth of AkHz (FWHM) and a Rabi 
frequency of 100/ci^z. The results from this modeling are 
shown in figure IHl along side the experimental data. De- 
spite the model not including any homogeneous dephas- 
ing the main features of the simulated process tomogra- 
phy (figure 13 closely matches those of the experimental 
data. In particular the rapid decay of the population 



component of the Bloch vector compared to the coher- 
ence components. Further simulations indicated that the 
decay rate of the population terms can be reduced by in- 
creasing the ratio of the Rabi frequency to the linewidth. 
A suitable criteria for when the application of the Bang 
Bang pulse sequence is useful for preserving arbitrary 
quantum states is when the decay rate of the popula- 
tion terms in the presence of the pulse sequence is slower 
than that of the coherence terms in the absence of the 
sequence. For the present case where T2 = O.865 the sim- 
ulation indicates that to meet this criteria it will be nec- 
essary to achieve a ratio of Rabi frequency to linewidth 
of VLrf /ooinh^^^^' There is limited capacity to increase 
the Rabi frequency of the driving field without the pos- 
sible excitation of off-resonant transitions. Therefore for 
the application of the Bang Bang pulse sequence to the 
mi = —1/2 ^ +3/2 transition to be useful it will be 
necessary to reduce the inhomogeneous broadening of the 
transition by a factor of ~10. 

The large change in magnetic field sensitivity as the 
Critical Point field is approached with no correspond- 
ing change in the inhomogeneous linewidth of the transi- 
tion suggests that the inhomogeneous broadening of the 
transition at the Critical Point field is not due to mag- 
netic interactions. The inhomogeneous broadening at the 
Critical Point field is probably due to strain within the 
crystal which is not intrinsic to the site and can be re- 
duced by refining standard crystal growing techniques. 
A reduction in strain broadening by over an order of 
magnitude has been achieved in analagous materials 20j 
through reducing the dopant concentration. Irrespec- 
tive of reaching the desired ratio of Rabi frequency to 
inhomogeneous linewidth, methods of designing pulse se- 
quences more robust to Rabi and detuning errors have 
been proposed [2ll[2^ . 

In conclusion the application of the Bang Bang pulse 
sequence demonstrates that dynamic decoupling tech- 
niques are applicable to correcting quantum errors on 
nuclear spin transitions of Pr^^'.Y2SiO^. This work re- 
alises very long hyperfine decoherence times, greater than 
30s, in a solid state optical A system suitable for quantum 
information processing applications. 
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